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ABSTRACT 

We report striking changes in the broadband spectrum of the compact jet of the black hole transient MAXI 
J 1836— 194 over state transitions during its discovery outburst in 201 1. A fading of the optical-infrared (IR) flux 
occurred as the source entered the hard-intermediate state, followed by a brightening as it returned to the hard state. 
The optical-lR spectrum was consistent with a power law from optically thin synchrotron emission, except when 
the X-ray spectrum was softest. By fitting the radio to optical spectra with a broken power law, we constrain the 
frequency and flux of the optically thick/thin break in the jet synchrotron spectrum. The break gradually shifted 
to higher frequencies as the source hardened at X-ray energies, from ~ 1 0 1 1 to ~4 x 10 13 Hz. The radiative jet 
luminosity integrated over the spectrum appeared to be greatest when the source entered the hard state during the 
outburst decay (although this is dependent on the high-energy cooling break, which is not seen directly), even 
though the radio flux was fading at the time. The physical process responsible for suppressing and reactivating the 
jet (neither of which are instantaneous but occur on timescales of weeks) is uncertain, but could arise from the 
varying inner accretion disk radius regulating the fraction of accreting matter that is channeled into the jet. This 
provides an unprecedented insight into the connection between inflow and outflow, and has implications for the 
conditions required for jets to be produced, and hence their launching process. 

Key words: accretion, accretion disks - black hole physics - ISM: jets and outflows - X-rays: binaries 
Online-only material: color figures 


1. INTRODUCTION 

Relativistic jets are ubiquitous features of black hole X-ray 
binaries (BHXBs). When a BHXB is accreting in a low/hard 
X-ray state (hereafter the hard state; for state definitions, see 
Belloni 2010), synchrotron radiation from compact jets is 
commonly detected as a flat or inverted radio spectrum and 
an excess of infrared (IR) emission (Corbel & Fender 2002; 
Fender 2006). In the canonical soft state, both the radio and IR 


* Based on observations collected at the European Southern Observatory, 
Chile, under ESO Program IDs 087.D-0914 and 089.D-0970. 


emission are quenched by at least a factor of a few hundred at 
radio frequencies (e.g., Corbel et al. 2000; Homan et al. 2005; 
Russell et al. 2011b), and jets may be prevented from being 
produced altogether in this state. 

The conditions required for launching jets and the mechanical 
launching process are fundamental unknowns in the field of 
accretion. The nature of the quenching/recovery of jets over 
state transitions provides essential constraints for jet models. 
However, to date, few studies have suggested how the jet might 
evolve over these transitions (Corbel et al. 2013a). During 
the hard-to-soft state transition, the IR/optical/UV emission 
from the jet appears to fade first at the start of the transition 
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(e.g., Homan et al. 2005; Russell et al. 2010; Yan & Yu 2012), 
while the flat radio spectrum persists as the X-rays soften 
then finally fades before bright optically thin radio flares are 
observed, often associated with relativistic ejecta (e.g., Fender 
et al. 2009; Miller- Jones et al. 2012). Similarly, on the reverse 
transition on outburst decline, the IR emission from the jet 
appears to return after the radio emission (e.g., Coriat et al. 
2009; Russell et al. 2010; Miller-Jones et al. 2012; Corbel et al. 
2013a). 

Phenomenologically, the detailed interplay between the well- 
documented pattern of changes in the inflow over the transition 
(Belloni 2010) and the outflow remains unclear. In the hard 
state, the flat radio spectrum breaks to an optically thin power 
law around the IR regime (e.g., Corbel & Fender 2002; Migliari 
et al. 2010). The break frequency varies between sources and 
even in time for the same source (Gandhi et al. 2011; Russell 
et al. 2013). It is unknown if or how the jet break varies during 
state transitions, but knowledge of this would provide essential 
constraints for models of how jets are produced. 

MAXI J1836— 194 was discovered in 2011 August as a new 
X-ray transient (Negoro et al. 2011). Within days it was clas- 
sified as a BHXB due to its X-ray temporal/spectral properties 
and its bright jet seen at radio-IR frequencies (Strohmayer & 
Smith 2011; Miller-Jones et al. 2011; Russell et al. 2011a; Reis 
et al. 2012). The source underwent a state transition from the 
hard state to the hard-intermediate state (HIMS; Belloni 2010), 
but instead of softening further it returned to the hard state and 
faded toward quiescence (Ferrigno et al. 2012). This type of 
“failed transition” outburst has now been reported in at least 
four BHXBs (e.g., Capitanio et al. 2009; Ferrigno et al. 2012; 
Soleri et al. 2013). 

Here we present the evolution of the jet spectrum of MAXI 
J 1836— 194 during its 2011 outburst using quasi-simultaneous 
multi-waveband radio, mm, mid-IR, near-IR, optical, UV, and 
X-ray data. 

2. BROADBAND OBSERVATIONS 

MAXI J 1836— 1 94 was observed with the Karl G. Jansky Very 
Large Array (VLA) from 2011 September 3 to December 3. 
We only consider data with quasi-simultaneous mm or mid/ 
near-IR data available. The radio observations were split across 
frequencies from 1 to 43 GHz. Due to the compact array 
configuration for most observations, we did not consider the 
lowest-frequency data (<4 GHz) as the target could not be 
easily resolved from a confusing source to the southwest. Data 
reduction was carried out according to standard procedures 
within the Common Astronomy Software Application (CASA). 
External gain calibration was performed using the primary and 
secondary calibrators, 3C286 and J1820— 2528, respectively. 
After phase-only self-calibration, the source flux densities were 
measured by fitting a point source in the image plane. Since no 
data were taken within three days of October 27, we estimated 
the radio spectrum on that date by fitting the smooth decay at 
each frequency with an exponential function, and interpolating 
for the appropriate MJD. 

Submillimeter Array observations were taken on 2011 
September 13 and 15 at frequencies of 256.5 and 266.8 GHz. 
We used 3C454.3 as a bandpass calibrator, Neptune to set the 
flux scale, and the phase calibrators 1911-201 and 1924-292. 
Data were reduced using standard procedures within CASA. 
MAXI J 1836— 194 was significantly detected on September 13 
(69.7 ± 6.9 mJy at 256.5 GHz; 66.4 ± 6.3 mJy at 266.8 GHz; the 
total on-source time was ~3.5 hr), but poor weather conditions 
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on September 15 prevented us from placing any meaningful 
constraint on the source brightness. 

Mid-IR observations of MAXI J1836— 194 were made 
with the Very Large Telescope (VLT) on four dates dur- 
ing its 2011 outburst. The VLT Imager and Spectrometer for 
the mid-infrared (VISIR) instrument on UT3 (Melipal) was 
used in small-field imaging mode. Observations were per- 
formed in 7347/7 (8. 2-9.0 /im), SIV (10.3-10.7 /rm), and .7/2.2 
(1 1.7-12.2 pm) filters and, on some dates, K hand (2. 0-2. 3 pm). 
Half of the approximately 1 hr observing time on each date was 
on source, due to the chop-nod mode. The data were reduced 
using the VISIR pipeline. Raw images were recombined and 
sensitivities were estimated based on standard star observations 
taken on the same night. For the K band, flux calibration was 
achieved using two K ~ 12.4 mag stars from the Two Micron 
All Sky Survey (Skrutskie et al. 2006), both lying ~6" from 
MAXI J 1836- 194. 

Optical images of MAXI J 1836— 194 were acquired on 
six dates spanning September 5-30 using the 2 m Faulkes 
Telescopes (FTs) North and South, located at Haleakala on 
Maui, USA, and Siding Spring, Australia, respectively, as part of 
an ongoing monitoring campaign of X-ray binaries (Lewis et al. 
2008). Images in Bessell B, V, R and Sloan Digital Sky Survey 
(SDSS) i' bands were reduced and calibrated using methods 
similar to those applied to other X-ray binaries observed by the 
FTs (e.g., Cadolle Bel et al. 2011; Russell et al. 2011b). Flux 
calibration in V and B bands was achieved using the Ultra Violet/ 
Optical Telescope (UVOT) calibration below. Landolt standard 
stars and a number of other calibrated fields in the R and i' 
bands were used to calibrate the field in these filters, adopting 
the transformation to SDSS i' band from R and 7 described in 
Jordi et al. (2006). 

In the X-ray band, Swift and the Rossi X-ray Timing Explorer 
(RXTE) monitored MAXI J1836— 194 every few days during 
its outburst (see Ferrigno et al. 2012). UVOT on board Swift 
observed in six optical/UV filters, uvotimsum was used to 
combine images in the same filter and uvotsource was used 
to perform aperture photometry. Magnitudes are based on the 
UVOT photometric system, which differs from the Bessell 
system by V — v < 0.04 and B — b < 0.04 for all reasonable 
spectral indices (Poole et al. 2008). We thus treat UVOT v, b 
and FT V , B magnitudes the same. Following the procedure 
outlined by Altamirano et al. (2008), we calculated colors and 
intensities using the 16 s time-resolution Proportional Counter 
Array Standard 2 mode data. For MAXI J1836— 194, we defined 
hard and soft color as the ratio of count rates in the 16-20 and 
6-16 keV bands, respectively, to that in the 2-6 keV band, 
and the intensity as that in 2-20 keV band. All values were 
averaged for each observation, and normalized to the Crab 
nebula on a per Proportional Counter Unit basis. 

3. RESULTS 

3.1. Light-curve Evolution 

The multiwavelength light curves are presented in Figure 1 . 
While the X-rays brightened, softened, then hardened and faded 
toward quiescence, the UV flux remained remarkably constant 
to within a factor of ~2 until November 14 (in quiescence 
the UV wl and w2 magnitudes are >3 mag fainter than this 
mean outburst level; Yang et al. 2012b). The optical B and V 
bands also remained approximately constant except during the 
intermediate states, when the source was slightly fainter. In the 
IR regime, this dip in flux was more pronounced, and the flux 
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Figure 1 . Upper panel: radio, IR, optical, and UV light curves of the 2011 
outburst of MAXI J 1 836—194. Where error bars are not visible they are smaller 
than the symbols. As well as VLT / VISIR mid-IR, FT optical, and Swift UVOT 
optical/UV data ( v band, b band, and the average flux of the three UV filters, w 1, 
m2, and w2) we also include i'- and Xs-band fluxes on August 30 (MJD 55803) 
from Rau et al. (201 1). Center panel: X-ray ( RXTE) light curve. Lower panel: 
X-ray color (two intensity ratios are given). The vertical dashed lines illustrate 
state transition dates (from Ferrigno et al. 2012). 

(A color version of this figure is available in the online journal.) 

was not constant after the transition. The amplitude of the flux 
drop/rise over the period of X-ray softening was a factor of > 10 
in the mid-IR filters, ~10 in K hand, ~2.5 in the optical i’ band, 
and <2 in B and the UV bands. The evolution is fairly smooth, 
with day-to-day rms flux variations of ~10%. 

While most BHXBs possess blue optical/IR spectra with 
flux densities peaking around the blue-optical/UV (e.g., Hynes 
2005), we find that the IR flux densities of MAXI J 1 836—194 are 
much brighter than the optical (the brightest we report is 57 mJy 
at 10.5 /iin). Very few mid-IR detections of transient BHXBs 
exist in the literature, and only GRO J0422+32, GX 339—4, and 
Cyg X-l (a high-mass BHXB) were detected at such bright flux 
densities (van Paradijs et al. 1994; Fender et al. 2000; Corbel & 
Fender 2002; Gandhi et al. 201 1). 

3.2. Broadband Spectra 

In Figure 2 we present broadband (radio to UV) spectra 
at six epochs in 2011 when radio and mm/IR observations 
were made within a few days of each other to constrain the 
jet synchrotron spectrum. Optical/IR data were dereddened 
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assuming an interstellar extinction of E(B — V) — 0.49 ± 0.03 
(Ay ^ 1.5 ± 0.1; this is inferred using broadband radio to 
X-ray spectral fitting which will be presented in a follow-up 
publication; T. D. Russell et al., in preparation) and adopting 
the extinction curves of Cardelli et al. (1989) and Weingartner 
& Draine (2001). 

The resulting broadband spectra are some of the most well- 
sampled in the literature for a BHXB. The radio and IR-optical 
spectral slopes can be approximated by power laws consistent 
with optically thick and optically thin synchrotron emission, 
respectively, confirming the dominance of the compact jet from 
radio to optical frequencies (see also Miller-Jones et al. 2011; 
Rau et al. 2011; Yang et al. 2012a). The radio power law 
appeared more inverted than most BHXBs (one similar source 
is XTE J1 118+480; Fender et al. 2001). The radio emission 
faded during and after the state transition, whereas the IR faded 
before brightening dramatically. This secondary maximum in 
the IR light curve is common to several BHXBs (e.g., Kalemci 
et al. 2005; Coriat et al. 2009). 

On September 16, when the X-ray spectrum was softest, the 
IR jet emission had faded considerably and a separate optical 
component was clearly visible, which is consistent with the 
irradiated accretion disk (T. D. Russell et al., in preparation). 
However, the /Gbaud detection, slightly brighter than the optical, 
implies that the jet was still contributing, and with the mid-IR 
upper limits, still allows us to constrain the optically thin spectral 
index (shown by dashed red lines in Figure 2). 

We fitted broken power laws to the broadband spectra to 
infer the flux and frequency of the spectral break (the results 
are given in Table 1). The fits used data up to B band on all 
dates except when the disk made a strong contribution to the 
optical emission (September 12 and 16). The UV data were 
excluded from all fits because the dereddened UV flux had 
much greater uncertainties due to the stronger influence of 
extinction and a significant but uncertain (and evolving) disk 
contribution. We estimate a maximum possible disk contribution 
of 30% in the B band on all dates except September 12 and 
16. The systematic uncertainty arising from this possible disk 
contribution is estimated and propagated into the fits, resulting 
in larger positive errors associated with the spectral index and 
inferred jet break frequencies (Table 1). Other processes such 
as synchrotron-self Comptonization cannot bias the power-law 
fits since their emission dominates at energies much higher than 
we fit here. 

A dramatic evolution of the jet spectrum is evident (Figure 2). 
The jet break resided at mm frequencies in the first four 
epochs, then shifted by more than two orders of magnitude 
to higher frequencies during the transition back to the hard 
state, being detected between our mid-IR and optical data on 
October 27. This shift caused an increase in the IR flux of the 
jet; the changing jet spectrum is responsible for the wavelength- 
dependent drop/rise in flux in the IR/optical light curves 
(Figure 1). The observed IR spectral index and the wavelength 
dependence of the drop/rise in flux over the transition are both 
inconsistent with the hot accretion flow model proposed by 
Veledina et al. (2013). 

3.3. Relation with X-Ray Hardness 

The X-ray hardness-intensity diagram (HID) from RXTE 
observations is presented in Figure 3. The overlaid circles 
represent the frequency of the jet break (left panel) and the 
total radiative jet luminosity integrated over the spectrum up to 
the B band assuming a source distance of 8 kpc (right panel). 
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Figure 2. Radio to UV (dereddened) spectra with broken power-law spectral fits (solid lines; see the text). The dashed lines indicate possible ranges of parameters. 


The horizontal bars denote the uncertainty on the spectral break frequency on each date. 




(A color version of this figure is 

available in the online journal.) 








Table 1 






Results of the Power-law Spectral Fitting, and Inferred Jet Breaks and Jet Luminosities 



Date 

MJD (UT) 

Optically Thick PL, 

Optically Thin PL, 

Break Frequency, 

FluXbreak 

Ljet 

(VLA) 

(VLA) 

a thick 

Cfihin 

Eb (Hz) 

(mjy) 

(ergs -1 ) 

2011 Sep 3 

55807.12 

0.70 ±0.10 

-0.60 ±0.10 

2.307 1 !® x 10 11 

426! 6 1 9 7 7 8 

5.21 42 j 2 ^ x 10 35 

2011 Sep 12 

55816.97 

0.20 ±0.01 

< -0.45 

>2.67 x 10 u 

>68 

(1.45-2.94) x 10 35 

2011 Sep 17 

55821.97 

0.17 ±0.02 

-0.69—0.35 

(0.410-5.08) x 10 11 

47-76 

(0.82-1.23) x 10 35 

2011 Sep 26 

55830.95 

0 6O +0 02 

U.OO—0 03 

—0 62 +0 08 
u -° — 0.03 

7.04f 4 * 1 'j S g x 10 11 

22U39 1 

4.3070.42 x 1035 

2011 Oct 12 

55846.01 

0.51 ±0.01 

—0 75 +0 06 
u - ' -0.05 

5.542°d 7 3 3 8 x 10 12 

1807 29 

7.632 1 1 - 5 o 4 9 x 10 35 

2011 Oct 27 

55861.00 

0.26 ±0.02 

—0 64 +012 
u -°^-0.04 

4.20+ 7 ( j 5 7 3 1 x 10 13 

262 1 5 7 

4.77T.Q 49 x 10 35 


Notes. For spectral indices, the convention F v oc v a is adopted. The error bars are 90% confidence limits. The integrated (5 x 10 9 Hz to 7 x 10 14 Hz) jet luminosities 
assume a source distance of 8 kpc. 


The frequency and flux of the jet break provide an estimate of the 
minimum radiative luminosity of the jet (e.g.. Corbel & Fender 
2002). We cannot constrain the frequency of the high-energy 
cooling break in the synchrotron spectrum, which occurs in the 
UV-X-ray regime, and we consider these estimates to represent 
lower limits on the radiative jet power. The total jet power will 
greatly exceed these values, as most of the power is locked in 
the kinetic flow of matter, and the jet radiative efficiency may 
be ~5% (e.g., Fender 2001). Using a composite Monte Carlo 
Spearman’s rank correlation test, we find tentative evidence for 
a positive correlation between the jet break frequency and X-ray 
hardness (as visually apparent in Figure 3). The Spearman’s rank 
coefficient is 0.79^o29> giving a significance level of 0 . 88^022 
for the correlation. 

4. DISCUSSION 

We have witnessed the gradual fading and recovery of a 

compact jet in the BHXB MAXI J 1 836—194, when it underwent 

transitions between the hard state and the HIMS. This suggests 


that the stability of the jet changes when the source moves away 
from the hard state, likely becoming less powerful. The results 
suggest the jet power is sensitive to the changing structure of the 
accretion flow between the canonical hard state and the HIMS. In 
the hard state, theoretically predicted scaling relations between 
inflow and outflow are observed empirically, at least in some 
BHXBs (e.g., Gallo et al. 2012; Corbel et al. 2013b; Russell 
et al. 2013 and references therein). These relations assume that 
a constant fraction of the matter being accreted is channeled 
into the jet, with total jet power, Pj et oc m jct oc m disk (where m 
is the mass accretion rate). In the HIMS the radio luminosity 
of MAXI J1836— 194 was not an accurate tracer of the total 
radiative jet luminosity because the latter decreased while the 
radio remained bright. The hard state scaling relations therefore 
no longer hold in the HIMS. The radiative jet luminosity (up to 
7 x 10 14 Hz) increased modestly with X-ray hardness (Figure 3). 
The ratio /«j et /ffidisk was therefore reduced in this state compared 
to the hard state. The jet power probably made up an increasing 
fraction of the total energy budget as the spectrum hardened 
and the outburst decayed (unseen changes in the frequency of 
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Figure 3. HID of the outburst, with marker size and color indicating the inferred jet break frequency (left; double circles indicate ranges inferred) and jet luminosity 
(right). The arrows represent temporal evolution. 

(A color version of this figure is available in the online journal.) 
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the high-energy synchrotron cooling break could affect this). 
The jet appeared less powerful in the HIMS but did not quench or 
recover immediately. If the jet switched off/on instantaneously, 
we would have observed a decrease/increase in its luminosity on 
short timescales — at least as long as the light crossing time over 
the photosphere or emitting region, and at most the adiabatic 
expansion timescale or the synchrotron cooling timescale of the 
propagating matter in the jet (the latter depends on the magnetic 
field strength, and is on the order of minutes for IR and days for 
radio frequencies, adopting B = 10 4 G; Miller-Joneset al. 2004; 
Gandhi et al. 2011). Instead, the jet power varied on timescales 
of weeks, the same as the evolving inner accretion flow. 

Compact jets are extinguished (followed by a bright, discrete 
ejection) either around the transition between the HIMS and the 
soft-intermediate state (the “jet line” in the HID) or before this, 
associated with a drop in the X-ray variability, and they tend to 
reactivate as the source hardens again in the HIMS (e.g.. Fender 
et al. 2009; Miller-Jones et al. 2012; Corbel et al. 2013a). The jet 
in MAXI J 1836— 194 was never extinguished completely, and 
the source never entered the soft-intermediate state, which is 
normally associated with the quenching of the compact jet and 
the launching of discrete ejections. The jet fading/recovering 
occurred in the HIMS, so the process was triggered by some 
activity of this state. Generally in the HIMS, the X-ray power 
law is softer, the disk more prominent, and the fractional rms 
variability is decreased, compared to the hard state (Belloni 
2010). The disk inner radius is smaller in the HIMS (Ferrigno 
et al. 2012), perhaps preventing as much matter from entering 
the jet nozzle, decreasing the mj et /wdisk ratio, in comparison 
with the hard state. Thermal disk winds are detected in the 
HIMS and most strongly in the soft state (Neilsen & Lee 2009; 
Ponti et al. 2012), so our results are consistent with the disk 
wind and jet accretion powers being inversely related; the wind 


may help quench the jet (see also Miller et al. 2012). Whatever 
the physical process is, these observations can provide vital 
constraints for models of the jet production mechanism. 

A further interesting constraint can be made from how the 
jet luminosity increase manifests itself over the HIMS to hard 
state transition. While the IR brightened considerably, the radio 
faded slightly and the jet break shifted to higher frequencies. 
These phenomena could be explained if the size scale of the 
jet base, the magnetic field strength or the particle pressure 
increased (Heinz & Sunyaev 2003). In this scenario, the IR 
flux and jet break frequency would increase while the radio 
emission would not be affected significantly, if the jet opening 
angle is constant. In addition, the gradual increase in the jet 
luminosity as the source hardened is consistent with the delay 
in the IR brightening after the transition is completed, seen in a 
number of BHXBs (e.g., Kalemci et al. 2005; Coriat et al. 2009; 
Russell et al. 2010; Corbel et al. 2013a). It is not clear why 
the radio spectrum evolved from one that was slightly inverted 
when the X-rays were softest (a ~ 0.2, where F v oc v“), to a 
very inverted spectrum as it hardened again ( a ~ 0.5). However, 
these results are consistent with the spectrum of compact jets 
becoming optically thin when they switch off (Fender et al. 
2004; Miller-Jones et al. 2012) and flattening when they switch 
back on (Corbel et al. 2013a, 2013b). 

For the first time, the detailed changes of the energetics of 
a relativistic jet have been witnessed over a state transition. 
In time these observations can place stringent constraints on 
time-dependent jet models, providing a step forward in our 
understanding of how jets are launched. These results also 
demonstrate that ground-based mid-IR monitoring of X-ray 
binaries can open up a new window into jet physics in these 
systems, especially when coordinated with other wavelengths. 
Very few submillimeter observations of BHXBs exist, but the 
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measured ~70 mJy submillimeter fluxes of MAXI J1836— 194 
and implied fluxes exceeding 0.1 Jy (Figure 2) suggest future 
outbursts of BHXBs will be easily detected by several current 
submillimeter facilities. 
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